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Qh! Abstract 

The contributions to the (g-2) of the electron from the corner-ladder graphs in 
sixth-order (three-loop) QED perturbation theory are evaluated in closed analytical 
^ ■ form. The results obtained are in excellent agreement with the most precise numerical 
evaluations already existing in the literature. Our results allows one to reduce the 
numerical uncertainty of the theoretical determination of the g-2 of the electron. 
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We have calculated in closed analytical form the contribution to the anomalous 
magnetic moment of the electron at sixth-order (three-loop) in QED perturbation theory 
from the so-called corner-ladder graphs shown in Fig.l. 

The result are, accounting for the mirror graphs (the infrared terms proportional 
to In A and In 2 A are omitted, they are listed in ref.[l]): 



a(graph 7) = -^(5) + ^(3) + ^tt 4 + 2tt 2 In 2 2 - 28a 4 - ? - In 4 2 - ^(3) 

55 2l 2743 2 2521 
+ 18 7T ln2 "l296 7r + W 

a(graph 9) =-C(5) - -tt 2 C(3) + — tt 4 + -n 2 In 2 2 - -a 4 - - In 4 2 
is f > 24 sv ) 72 si ) 2160 g 3 18 

107 25 2n n 3 2 43 

C 3 H tt 2 In 2 H tt 2 H , 

12 sv ; 18 16 36 

a(graph 17) =| C (5) " jj- 2 C(3) - |U 4 - H^ 2 In 2 2 + ^a 4 + g In 4 2 

+ ^ C(3 )-^ln2 + ^ 2 -^, 
24 SV ; 18 648 27 ' 

a(graph 19) = ^C(5) - |tt 2 C(3) - -^tt 4 + ^tt 2 In 2 2 - ^a 4 - ^ In 4 2 
' 24 w 8 432 18 3 18 

635 83 2l n 4777 2 1835 

C(3) + — 7r 2 ln2 tt 2 H , 

72 sv ' 18 2592 864 ' 

, 215 ^ , 95 2>/ , 41 4 137 2l 2 160 20, 4 

a(graph 27) = C(5) + 7^ C(3) + — tt 4 - — tt 2 In 2 2 + — a 4 + — In 4 2 

69 _ 101 2n „ 2401 2 3017 

H C 3 tt 2 In 2 H tt 2 , 

4 sw 18 2592 864 ' 

a(total) =^C(5) + ^ 2 C(3) + ^tt 4 - 3tt 2 In 2 2 - 12a 4 - V 2 - i? C(3 ) 

31 2, 655 2 481 
tt 2 In 2 -| 7T 2 + 



(1) 



(2) 



(3) 



(4) 



(5) 



(6) 



6 216 288 

Here and in the following Q{p) is the Riemann ^-function of argument p, £(p) = 

— , (we recall that C(2) = tt 2 /6, C(3) = 1.202 056 903..., C(4) = tt 4 /90, C(5) = 



riP 

n — 1 



1.036 927 755...), and a 4 = V 7 = 0.517 479 061 



, 2 n n 4 

n=l 



The contributions of these graphs had been previously calculated only by (approxi- 
mate) numerical methods. The numerical values of eqs.(l)-(6) turn out to be in excellent 
agreement with the values shown in ref.[l] and [2] (see Table I). 
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If C3 is the coefficient of the ( — J term in the perturbative expansion of the electron 
anomaly in QED 

aMm ^)^)\^)\^)\..., (7) 

by using eq.(6) we have 

c 3 = 1.176 19(21) , (8) 

to be compared with the value [3] 

c 3 = 1.176 13(42) ; (9) 

the error of eq.(8) is due to numerical uncertainty of the value of the contribution of 
the remaining group of five graphs still not known in analytical form [4]. 
Using the experimental determination of a [5], 



oT^exp) = 137.035 997 9(32) , (10) 



the values of C2 and C4 [6] 



197 1 3 

° 2 = 144 + 2 C(2) + 4 C(3) " 3C(2) ln2 = "°- 328 478 965 - 

c 4 = -1.434(138) (11) 

and accounting for the small vacuum polarization contributions of muon, tau and hadron 
loops as well as the electroweak contribution [6] one finds 

a e (theory) = 1 159 652 141.4(27.1)(2.6)(4.1) x 10" 12 (12) 

where the first error comes from the error of the determination of a (10), the second 
from the value of C3 and the third from the error of C4. The error due to the three- loop 
coefficient C3 is now less than the error due to C4. 

If one assumes the validity of QED, using the experimental value of the electron 
anomaly [7] 

a e (exp) = 1 159 652 188.4(4.3) x 10" 12 , (13) 
one can estimate a from eq.(7); one finds 

« _1 (a e ) = 137.035 992 34(51) (31) (48) , (14) 

where the errors come respectively, from eq.(13), from C3 and from C4. Summing up the 
errors of eq.(14) in quadrature one finds 

a-\a e ) = 137.035 992 34(77) (15) 
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which is the most accurate determination of a available at present. 

We describe now briefly the techniques used for deriving eqs.(l)-(6). All the vertex 
graphs of Fig.l can be obtained from the self- mass-like graph of Fig. 2 by inserting the 
external photon line in all possible ways along the electron line. Once that the limit 
A — >■ is taken, the contribution from each of graphs of Fig.l is expressed as a sum 
of some hundreds of terms all of the form N/D, where iV is a scalar polynomial in p 
and the internal integration momenta, and D = D™ 1 D^ 2 ... -Dg 8 (n^ = 0, 1, 2), the 
being the denominators of the various propagators of the graph of Fig. 2. The simplest 
term, with N = 1 and all rii = 1, was calculated in analytical form in a precedent paper 



The contributions of graphs 19 and 27 turn out to be finite; on the contrary, 
graphs 7, 9 and 17 need renormalization, and as consequence their contributions become 
infrared divergent; in this case appropriate ultraviolet counterterms are subtracted, as 
well as infrared ones which render the entire contributions finite. 

Let us consider for instance the term defined as 



with N = (k 1 - k 2 ) 2 and 

D={k\- ie){kl + 1 - ie)((fci - k 2 f + 1 - ie)((p - h) 2 + 1 - ie)((p - k 2 ) 2 - ie)x 
x(q 2 + 1 - ie)((q - h) 2 + 1 - ie)((q - k 2 ) 2 - ie) 



(the electron mass is set equal to 1). In order to evaluate eq.(16) we apply the techniques 
used for calculating the analytical contribution to the g-2 of the muon of the "light- 
light" vertex graphs [9]. The self-mass-like graph of Fig. 2 is depicted so as to show 
the topological similarity with the correspondent self- mass-like graph obtained from 
the "light-light" graphs. Following [9], we write a dispersion relation in (ki — k 2 ) 2 at 
constant k\ = I, k\ = m, for the inserted g-loop vertex part 



[8]. 




(16) 



/ 



(-i)d 4 q 



(q 2 + l- ie)((q - h) 2 + 1 - ie)((q - k 2 ) 2 - ie) 
o f°° db 1 




b + (h - k 2 ) 2 R(b,-l,-m) 



where 



In W ~ In b ^ + 1 + ^ + ^ - 1 + ^ + ^ - ^ R ( b ' ~ 1 ' 
n n b(b + I + m) + (b - I + m) - (b - 1) R(b, -/, -m) 

and R(x, y, z) is the usual two-body phase space square root 



(18) 



R(x,y,z) = a/(x — y — z) 2 — Ayz = \/ x 2 + y 2 + z 2 — 2xy — 2yz — 2xz . 



(19) 
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Now we introduce hyperspherical variables for both the k\ and the k 2 loops 

d 4 ki = -I dl dQ^ki) , 

1 „ ( 20 ) 

d A k2 = -to dm dfl 4(^2) , 

and perform the hyperspherical angular integrals by expanding the denominators in 
Gegenbauer polynomials [10]. In the mass-shell limit p 2 = —1 we obtain 



/ 



dQ^dQ^k) _ (2tt 2 ) 2 iny± ^ 



((p - &1) 2 + l)(p - A; 2 ) 2 (6 + (fci - A; 2 ) 2 ) /m 
where 

1 



r± = i + 



8m/ 



v / /(/ + 4) -/ [(m+ 1) ± (to- 1)] [6 + Z + to- R(b,-l,-m) 



The analytical continuation for timelike values of p 2 of the radial integral over m requires 
also a deformation of the contour [10] of the integration in order to avoid the singularity 
at m = —1. The contour is split in two parts, one where m > — 1 and another where 
— 1 < m < 0. Y- must be used when integrating over m in the first zone and and Y + in 
the second one. At this point the integral (16) reads 

I = r*(rjn> r _ g ; lnWlnY _ 



\nWlnY + ) . (22) 



x m+1 J 1 R(b,—l,—m) 

r° dm r°° db 

7_i to + 1 J 1 R(b, -I, -to) 



Eq.(22) contains two square roots: R(b, —I, —to) and a//(/ + 4). We choose to integrate 
eq.(22) in b first, then in m and at last in /; we used this order of integration for all the 
other integrals over 6, m and / needed in the calculation of eqs. (l)-(6) 1 . If we write 



f°° dl 



the b and to integrals are performed by differentiating (repeatedly, when needed) with 
respect to / (see ref.[8]), then integrating the resulting (simpler) expression, and re- 
covering f(l) by a quadrature. Note that the derivatives with respect to /, to and 



1 Note the simplification with respect to the analogous "light-light" graphs calcula- 
tion; in eq.(23) of ref.[9] four square roots appear, and as consequence the decomposition 
of the hyperspherical logarithm is needed and different orders of integration must be 
chosen for the different pieces. 
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b of the functions \nW and lnY"-j- contain the same polynomial in the denominator, 
(b + m) 2 + l(b — l)(m + 1), even if the two functions have different origin; this fact 
simplifies considerably the calculations. At last, f(l) is expressed as a sum of known 
tetralogarithms of the variable x = (I + 2 — a//(/ + 4))/2 : 

/(/) = - 9Li 4 (x) - 15LU(-x) + \s 2 , 2 {x 2 ) - 2S 2 , 2 {-x) + 2T 2 , 2 (x) - U 2 (x) 

+ lux (5Li 3 (x) + 8Li 3 (-x) - 2S h2 (x) - 2U 2 (x) ln(l - x)) (23) 
- In 2 x (u 2 (x) + hi 2 (-x) + X - ln 2 (l - x)\ + 2((2)U 2 (-x) , 



where 



and 



( — l)n+p-l r 1 fa 

S n , p (x) = [ {n \ )lpl J o j ^(t) h?(l - xt) , 
Li n+ i(x) = S n) i(x) , 

T 2 , 2 (x)= / ^ln(t + l)Li 2 (t) . 
Jo t 



The last definite integral over / is performed using known tables of integrals containing 
poly logarithms; one finds 

/ = yC(2)C(3) - f C(5) • (24) 

The processing of the other terms with different iV and D is similar to that described 
above. When the numerator iV contains powers of q we use the decomposition described 
in ref.[9]; as this fact causes a severe blowing-up of the size of the calculations, some 
algebraic manipulations are used in order to reduce the maximum power to q 3 . In 
the course of the calculations, some non-standard functions of polylogarithmic degree 4 
appear, as an example 



gx(x) = [ dt\n(l-t + t 
Jo 



2) Li 2 (t) 



the expressions of the functions gi(x) by means of polylogarithmic functions S njP (x) are 
not known. The last integration over / of the functions gi(x), after some integration by 
parts, give a group of definite integrals containing ln(l — t + 1 2 ) which we found to be 
the same ones appeared in the calculation of the contribution to the g-2 of the electron 



2 The expression of the function T 2j2 (x) in terms of polylogarithmic functions S njP (x) 
is not known; nevertheless, the properties of T 2y2 (x) are similar to those of S niP (x), so 
that we consider this function as a "standard" polylogarithm. 
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of another set of vertex graphs [11], and which were calculated in analytical form in 
ref.[12]. 

The final infrared-cutoff dependent results (l)-(6) are obtained using the analytical 
expressions of the infrared counterterms shown in ref.[13]. 

As a check, for each graph of Fig.l we have modified the internal routing of the 
momentum A and we have extracted the contribution to the (g-2), obtaining new ex- 
pressions, whose subsequent analytical evaluations reproduce correctly eqs.(l)-(6). 

The algebra of the whole calculation was processed relying on the algebraic manip- 
ulation program ASHMEDAI [14] . All calculations were done on a number of VAXsta- 
tion 4000/60 and 4000/90. The processing of the whole contribution of a graph required 
about 40h of continued CPU time and 150 Mbytes of disk space for storing of interme- 
diate expressions on a VAXstation 4000/90, the peak length of the algebraic expressions 
being about half million of terms. Due to the (relative) high speed of machines used, 
we devoted our efforts to design and carry out the maximum quantity of consistency 
and cross checks, rather than to speed up the execution time or to minimize the disk 
space requirements. 

Work at the contribution to the electron anomaly of the remaining graphs still not 
known in analytical form is now in progress. 
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Figure captions 

Fig.l: The corner-ladder graphs (mirror graphs are omitted). The numbering follows 
ref.[l]. 

Fig. 2: The self- mass-like scalar graph with the cut of the dispersion relation. 
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Table I 



Graph 


r r-i 1 

ref.[l] 


ref. [2] 


Our 


7 


-2.6707(19) 


-2.670 546(30) 


-2.670 554 745 651... 


9 


0.6189(64) 


0.617 727(121) 


0.617 711 782 178... 


17 


0.6097(34) 


0.607 660(240) 


0.607 752 806 792... 


19 


-0.3182(72) 


-0.334 698(11) 


-0.334 695 103 723... 


27 


1.8572(86) 


1.861 992(240) 


1.861 907 872 591... 


Total 


0.0893(60) 


0.082 065(362) 


0.082 122 612 187... 



Table I: Comparison of our values of the contributions of the graphs shown in Fig.l 
with the values of ref.[l] and ref. [2]. 



10 



This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9410248v2 



